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It is usually argued that the presence of gapless quasiparticle excitations at the nodes of the d- 
wave superconducting gap should strongly decohere the quantum states of a d-wave qubit, making 
quantum effects practically unobservable. Using a self-consistent linear response non-equilibrium 
quasiclassical formalism, we show that this is not necessarily true. We find quasiparticle conductance 
of a d-wave grain boundary junction to be strongly phase dependent. Midgap states as well as 
nodal quasiparticles contribute to the conductance and therefore decoherence. Quantum behavior 
is estimated to be detectable in a qubit containing a d-wave junction with appropriate parameters. 



Among numerous qubit implementations, supercon- 
ducting ones enjoy long decoherence times because of 
their gapped electronic excitation spectrum. This fact 
has recently been confirmed by several striking experi- 
ments [3 . The key constituents in all of those are Joseph- 
son tunnel junctions. One advantage of the tunnel junc- 
tions is the exponential dependence of their quasiparticle 
resistance R on temperature T (i.e. R ~ e^/^ where A 
is the superconducting gap. Herein kB=fi=^)- The elec- 
tronic decoherence is therefore exponentially suppressed 
at low T. Similar behavior also exists in superconduct- 
ing point contacts, with the energy of Andreev levels 
eo{(j)) = A cos 0/2, replacing A in the exponent 0, 13- 
Deviation from the exponential dependence is expected 
at low temperatures 3] . 

Despite their naturally degenerate ground states 
desirable for quantum computation, d-wave qubits 0, 
are controversial because their quasiparticle spectrum is 
gapless at the nodes of the order parameter. Moreover, 
experimentally, the normal resistance extracted from I- 
V characteristics of d-wave grain boundary junctions is 
found to be very small 7, 8] (with RC ^ Ips 8]); much 
smaller than required to observe quantum effects. How- 
ever, the resistance is measured in the running state 
of the junctions. The Doppler shift due to large su- 
perconducting current in such a state will populate the 
nodes, enhancing the conductance G (= 1/R). More- 
over, time variation of the phase difference across the 
junction would effectively phase average G. As we shall 
see, midgap states (MGS) can significantly contribute to 
such an averaged G, except at very low T. 

The first step is therefore to calculate G for a d-wave 
grain boundary junction. Most existing methods can 
study ac properties of a Josephson junction biased with a 
constant voltage (see e.g. 0,113] )• This however implies a 
constant variation of phase difference </), which is not the 
case in qubits. In Ref.|3, a self-consistent non-equilibrium 
quasiclassical technique was developed to calculate linear 
response of a Josephson junction to an ac voltage with 
arbitrary frequency. The method was successfully ap- 
plied to the case of a superconducting point contact. In 
high Tc superconductors, the relatively large Tc/Ep {Ep 
is the Fermi energy), makes the quasiclassical approxi- 



mation only marginally applicable. It nevertheless has 
proven successful in calculatirig equilibrium properties of 
d-wave superconductors [ill Il3. fl3j |. Here, we employ 
the theory of Ref. to calculate G, and therefore deco- 
herence, in a d-wave grain boundary junction. 

Let us now briefly describe the technique. More details 
are provided in Ref. 0. We calculate quasiclassical re- 
tarded, advanced, and Keldysh Green's functions 0, 
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which are functions of the Fermi velocity wp, quasi- 
particle energy e, position r, and time t. Here 
f''^{vp, e; r, t) = f^{~\-p, — e; r, t)* , etc. In equilibrium, 
the retarded and advanced Green's functions can be writ- 
ten in terms of Riccati amplitudes [l3| Oq ^'^d ^ 
very similar to the conventional method for the Matsub- 
ara Green's functions [l^ll3| : 
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where a = R, A for retarded and advanced functions re- 
spectively, and s" = + {—) for a — R (A). The subscript 
"0" denotes equilibrium quantities. The amplitudes sat- 
isfy Riccati-type equations 
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where e" = e + is°'r], with e and 77 being the real 
and imaginary parts of the quasiparticle energy respec- 
tively. In d-wave superconductors, 77 results from both 
inelastic and impurity scattering processes 0|. The 
boundary conditions are the bulk solutions of o,^ = 
Ao/(- ze" + s"17") a nd = AS/(-ie" + s"!)"), where 
r2"=v/|Aop - (e")2. To calculate and {b§ and 
Qq), we integrate (jJJ in the direction of vp {~vp), start- 
ing from the boundary conditions at x— — 00 (-t-00). 

Non-equilibrium behavior emerges when an ac voltage 
V{t) = Vocoswt is apphed across the junction. Choos- 
ing a gauge in which the vector potential A = 0, we 
consider a scalar potential <& — ±(Vb/2) coswi, with the 
+ (— ) sign on the left (right) side of the junction. To 
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ensure applicability of the linear response formalism, it 
is necessary that eVb <C A,w [l5| . 

We define the linear response Green's functions 6g°' = 
5" ~ ?o J related to the linear response amplitudes da" = 
a" - and Sb" = b" - through (a = R, A) 

6a"b?; 

Sg" = -2s"- 
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where Uq^. = flg (e ± uj/2), etc. We also introduce the 
anomalous Green's function Sg'^ (with the same matrix 
form as 5^) by (5g^ = Sg^ {T+ 5g^T- - Sg^T+ , 

where JF± = tanh[(e ± ti>/2)/2T]. Correspondingly, we 
introduce anomalous functions Sa'^ and (56^ which are 
related to the Green's functions through 
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The differential equations describing Sa" are Vi?-V(5a" = 
A"Sa" + B", with bulk boundary conditions Sa" = 
-B"/A", where 
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The corresponding equations for Sb" can be obtained by 
applying a f-operation to both sides of the above equa- 
tions. Integrations for Sa^, Sb"^, and Sa-^ [Sb^, 5a^ , and 
Sb^) are taken in the direction of (— v^), along the 
quasiclassical trajectories. 

The equilibrium order parameter, Aq, is calculated self 
consistently using numerical iteration, while the linear 
response (5A = A — Aq is given by 



^A(w) = -2e(Ao/w)$(w). 



(7) 



(5A satisfies the self-consistency relation up to an error of 
0(jc/jc,buik), where jc O'cbuik) is the Josephson (bulk) 
critical current density. In d-wave grain boundary junc- 
tions, jc/jc,buik is smallest for a 0°— 45° junction. Appli- 
cation of O, instead of common self-consistent iterative 
methods 16], is a clear advantage of the present tech- 
nique. 

We consider an asymmetric 0°— 45° misoriented d-wave 
grain boundary junction, with perfect transparency and 
no roughness. To clearly see the effect of the nodes and 
bound states, it is useful to define an angle resolved con- 
ductance Ge = RefJje/VbjS', where S is the the area of 
the junction and 
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FIG. 1: (a) Angle resolved conductance, normalized to Go = 
NfVfS , for (j) — Q (dashed line) and </> = 7r/2 (solid line) 
at T = O.OlTc. (b) Quasiparticle trajectories crossing the 
grain boundary. Trajectories along a and /3 directions see 
sign change of the order parameter. Other parameters are: 
Lo = 0.02rc and ?? = O.OSTc 



is the contribution to the ac current density from a tra- 
jectory in 9 direction (cf. Fig. QJd). Here % is the third 
Pauli matrix, n is a unit vector in the direction of the 
trajectory (n • x = cos 61), Np is the density of states at 
the Fermi surface, and tc is an energy cutoff, which in our 
numerical calculations is taken to be lOTc. (The results 
however are insensitive to the exact value of ec ) The 
total conductance is given by G = ^ {dO 1211)08 cos 6*. 

Figure shows Ge for two values of phase difference 
<f) across the junction. At = 0, G is maximum a,t 9 = 
-22.5° and 67.5° (paths a and (3 in Fig.UlD). These are 
the directions along which the quasiparticles see different 
signs of the d-wave order parameter on either side of the 
boundary. Zero energy Andreev bound states formed in 
these directions are responsible for the large G ||3| . Notice 
the asymmetry of Go with respect to = at (/) = 0. 
This asymmetry leads to an ac current parallel to the 
boundary The current vanishes at (j) = 7r/2 and is 
not dissipative (it is perpendicular to the electric field), 
therefore not central to our discussion. 

A nonzero phase difference splits the bound states, 
moving them away from zero energy. Maximum split- 
ting occurs aX (f) — n/2. Nodal directions dominate the 
quasiparticle current at this phase difference (cf. Fig.^). 
Notice orders of magnitude difference in the conductance 
maxima between the two phase differences. 

Figure displays (/)-dependence of the total G at two 
temperatures. Sharp peaks at = and tt, 2-3 orders 
of magnitude larger than the minimum (at (p — 7''/2), 
result from the MGS. The conductance spikes are much 
narrower at T = O.OlTc than at O.lTc and become ex- 
tremely sharp at even lower T. As a result, the infiuence 
of the MGS is small at low T (see below) . Figure 
shows the T-dependence of G at </> = and it/2. While 
the two curves coincide at high T, they diverge as T is 
lowered. The conductance at (/> = 7r/2 is suppressed at 
low T . This is because the gapless quasiparticles on one 
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FIG. 2: <j> (a) and T (b) dependence of G for the same set of piQ 3. Dependence of G on lj (a) and rj (b) at T = IQ-^n. 



parameters as in Fig. Q 



77 in (a) and u> in (b) are equal to O.OlTc 



side of the junction see a gapped order parameter on the 
other side As T — > 0, however, it starts to saturate 
to a value almost proportional to 77 (see Fig.lJla). Similar 
behavior was also found for the case of superconducting 
point contacts Q. 

To reduce decoherence, the qubit should be operated 
at low T. Figure O shows the u and 77 dependence of G at 
T = IO^^'Tc. While strong w-dependence of G exists at 
(/i = 0, it is almost frequency independent when = 7r/2 
(Fig. The //-dependence of G at (fi — tt/2, on the 

other hand, is close to linear (Fig. |31d). The positive 
slope can be attributed to the broadening of the MGS 
with 77, which enhances their overlap near zero energy, 
increasing their contribution to G. 

We now try to calculate the decoherence time in a d- 
wave qubit, assuming a double well potential U{(l)) with 
minima at ztcj)^ « ±7r/2. (A practical example of such 
a system is given in Ref. |6|.) We first write down a 
Hamiltonian which, in classical regime, reproduces the 
above calculated G. At low frequencies, G shows slow uj- 
dependence for almost all (j) (except for (/) = 0, tt, which 
do not contribute to decoherence — see below). One can 
therefore assume our system to be coupled to an Ohmic 
heat bath 0. To get a (/)-dependent G, however, it 
is necessary to consider nonlinear coupling to the bath 
[H 1131 [see Eq. below]. We therefore write the 

Hamiltonian as 

+ UW ~ F{<jy)X, (9) 

where Q = 2eP^ {P^ momentum conjugate to (p) is the 
charge stored on the junction, and X is a heat bath oper- 
ator. To find the relation between the nonlinear coupling 
function F and G, we first study the classical behavior 
of the above Hamiltonian. The equation of motion is 

(G/4e2)0 = ~d^U{cp) + d^F{c^)X, (10) 

Assuming small back-action of the system on the heat 
bath, one can use linear response theory (20,] to write 

X{t) = xW(t) + / dt'D{t - t')F[ci,{t% (11) 



where X'"^ is the unperturbed bath operator, D{t — 
t') = (i[X(o)(t),X(o)(t')]-)6l(t-i') is the bath's retarded 
Green's function, and (...) denotes equilibrium statistical 
averaging. For an Ohmic heat bath D{t) = —a5'{t) [in 
Fourier space D[u)) = zaw], with a being the dissipation 
coefficient and 6'{t) defined by / dt'6'{t - t')z(t') = z{t). 

Eq. ((Tn|l therefore leads to a stochastic equation that 
resembles the classical Langevin equation 

iC/Ae^)4> + a{d^Ff,j> + d^U = ^^, (12) 

where £,4> = id^F)X'^^'> (t) is the fluctuation force [2l| . 
The second term in (|12ll gives dissipation and there- 
fore corresponds to the 0-dependent G G{(j>) = 
Ae'^a {d^F)^ . As a result 

F{<ty)^-^ f dc^'^GW). (13) 
2eVa Jo 

Notice the nonlocal dependence of F on G. The spike 
in G at = (cf. Fig. Et) becomes extremely sharp 
as T^O, with little contribution to F{ip). This, although 
easily justified numerically, can be understood by making 
an analogy to the case of a superconducting point contact 
0: G(0) - (A/T)e-'^l'^l/^ near the spike {(j) measured 
from the center of the spike) , leading to a contribution to 
F(0) proportional to -y/T/A, which is negligible as T^O. 
Numerical integration shows an even smaller effect. 

In the quantum regime, one can truncate the Hilbert 
space to the left and right degenerate states, 
of the double-well potential, with {L\F{(j))\R) ~ and 
{R\F{(I))\R) ~ -(L|F(0)|L) ~ F(0o). Equation © then 
gives the effective two-state Hamiltonian of the system 
(at the degeneracy), coupled to the heat bath, as 

H,s = -^a, - Fi<j)o)Xa„ (14) 

where SE is the energy splitting between the two low- 
est energy states of the system. The dephasing rate is 
proportional to the spectrum of heat bath fluctuations, 
S{uj) = Im[Z)((jj)] coth(ti;/2T), taken at the resonance fre- 
quency of the two-level system [see, e.g., Eq. (3.11) in 
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Then, the dephasing time due to coupHng to 
the Ohmic heat bath is obtained from 

6E 

T^^ =aF{MHEcoi\,—. (15) 

Numerical calculation at T = lO^^Tc gives j^^"^ \fGd(t) w 
0.38v^, for 77 = O.OlTc. Substituting into ^ and ^ 
and restoring h, we find w 0.023RqGo5E /h, where 
^Q=^/(2e)^«6.45 kl7 is the quantum resistance. 

For quantitative estimation of t^, we need to know 
the value of Gq. We extract Go from the Josephson 
critical current density whose value is available from ex- 
periment 0: jc ~ 10^—10"' A/cm^. Our calculations 
at T = 0.05Tc (close to 4.2K where most experiments 
are performed) show a critical current Ic ~ 0.08GoTc/e, 
almost independent of rj. For a submicron junction of 
area S ^ 0.01 /xm^, we obtain (taking Tc ~ 100 K) 
Go - 10-^-10-3 13. Assuming SE/h - 1 GHz, we 
find ^ 1—100 ns (qubit quality factor Q ^ 1—100). 
Smaller 77 will result in a larger r^. It is therefore de- 
sirable to use materials with low disorder and junctions 
with small roughness. Depending on the parameters, it is 
possible to observe signatures of quantum behavior (e.g. 
coherent tunneling) , although the decoherence time may 
not be long enough for quantum computation. 

Realistic junctions, in general, suffer from finite reflec- 
tivity, roughness, and faceting [23. The effect of rough- 
ness, to some extent, is similar to that of 77; it broadens 
the MGS, increasing their contribution to G as they split 
■ Imperfect transparency also affects the MGS in a 
nontrivial way, influencing G. Presence of a subdominant 
order parameter at the junction, although yet unjustified 
experimentally [2^ . can enhance the quasiparticlc resis- 
tance by opening a gap at the nodes. A small bulk size 
can also have a similar effect by quantizing momentum 
along the nodal directions. Most of the above mentioned 
effects can be studied within the framework of the present 
formalism and are the subject of further investigations. 

Other sources, such as fluctuations of the external 
fields, coupling of bulk nodal quasiparticles to the electro- 
magnetic field produced by the qubit (due to e.g. spon- 
taneous currents), coupling to nuclear spins or param- 
agnetic impurities, background charge fluctuations, etc. 
will also contribute to the decoherence. 
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